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ABSTRACT 

We report the serendipitous discovery of transient, point-like high energy 7-ray emission coincident with the position 
of the suspected BL Lac object 5C 3.178. The source was detected using the Fermi Large Area Telescope (LAT) at a 
significance level of 8a during a 200 day period which began on November 17th 2011 (MJD 55882). The observed 
7-ray emission is characterised by a moderate 0.2-300 GeV flux -Fb. 2-300 Gcv = (8.22 ± 2.04) X 10"^ ph cm s and 
a hard power law spectrum with spectral index F — 1.76 ± 0.09. These properties appear consistent with the known 
sub-population of TeV 7-ray-emitting BL Lac objects, implying that the source may be detectable using atmospheric 
Cherenkov telescope arrays. Moreover, the temporal variability of the source suggests that real-time searches of the 
Fermi-LAT all-sky dataset for weak emission on ~ 200 day timescales may represent a rewarding strategy. 
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1. Introduction 



Oh, 

Q Active Galactic Nuclei (AGNs) are compact, energetic as- 
5— I jtrophysical phenomena associated with the cores of many 
galaxies, which are paradigmatically assumed to be pow- 
d ered by disk accretion onto a supermassive (M > 10^ A/©) 
black hole. A sizeable fraction of AGNs are also launch sites 
for relativistic jets of plasma, which emit broad-band radi- 
^ ation with non-thermal energ y spectra that ext end from 
QQ radio to 7-ray frequencies fe.g. lUrrv fc Padovanil [19951 ). 
00 ■ BL Lacertae objects (hereafter BL Lacs) are a subset 
CO pf AGNs that have jet axes that are oriented close to the 
CO observer's line of sight. The observed spectral energy dis- 
tributions (SEDs) of BL Lacs are dominated by relativisti- 
cally beamed emission from charged particles entrained in 
CO the jets and exhibit a characteristic two-component struc- 
y—i ture. The individual components are visible as low and high 
^ energy peaks in a vFi, representation of the SED and are 
. ^ typically identified with distinct radiative processes occur- 
ring within the jets. The low energy component is usually 
\—{ attributed to synchrotron emission while inverse- Compton 
. ^ .Up-scattering of this synchrotron emission or other ambi- 
ent soft radiation fields is assumed to dominate at higher 
energies. 

BL Lacs are an established class of variable high energy 
(HE; 10-3 < E^/GeY < 50) and very high energy (VHE; 
10~3 < E^/TeV < 10^) 7-ray sources. Indeed, the sec- 
ond Fermi-LAT (2FGL) source catalogue includes 435 con- 
firmed BL Lac ob jects, of which 36 h ave also been detected 
at TeV energies (iNolan et al.ll2012l) . In a recent review, 
iSentiirklTal] (|20ll) " examined the Fermi-LAT light curves 
of 12 established VHE 7-ray-emitting BL Lacs, which re- 
veal significant flux variability on timescales ranging from 
days to months. Furthermore, a population study of 395 BL 
Lac objects listed in the second Fermi-LAT AGN catalogue 



(2LAC) also revealed clear evidence for HE variability, al- 
though it is a pparently less prevalent in sources with hard 
7-ray spectra (| Acker mann et al"1l201lh . 

5C 3.178 is a point-like radio source located to the 
southeast of M3 1 at (a.T2nnn , ^.lannn) = (00M7'"55^220, = 
-t-39°48'57".60) ()Maol 120111 ). The source was initially 
discover ed using the Cambridge one-mile radio telescope 
(jPoolevI [l969.) and has been provisionally identified as a 
BL Lac object based on the absence of detectable emis- 
sion f eatures in its optical s pectrum (jBurbidge fc HewittI 
Il987t iDiorgovski et al.l Il995l ) . Multi-wave length counter- 
parts have been i dentified at near infra-re d (Skrut skie et al.l 
l2006l: lMaoll2011h. optical (IDiorgovski et ah 1995,). u ltravi- 



olet (iFittingoff et al.l [2OO9O and X-rav (|Voges et all 119991 : 
ISaxton et al.l 20081 ) frequencies. A source redshift of z = 
0.2517 ± 0.001 has been established on the basis of op- 
tical absorption featur es imprinted by the host galaxy 
(IDiorgovski et al.lll995[ ). 

The Large Area Telescope is the principal scientific in- 
strument on the Fermi Gamma Ray Space Telescope space- 
craft. It is a pair conversion telescope with a ~ 1 m effective 
area, which is sensitive to photons with energies between 20 
MeV and 300 GeV. The field of view of the Fermi-LAT en- 
compasses 2.4 steradians, and the instrument typically op- 
erates in a continuous survey mode, providing all-sky cov- 
erage every two orbits. The operation and performance of 
the F ermi-LAT are described in detail by (jAtwood et al.l 
[2OO9I) . 

This article describes the discovery and analysis of a 
point-like 7-ray signal coincident with the position of 5C 
3.178, which was serendipitously identified while search- 
ing unsuccessfully for HE emission associated with a 200 
day flare of the nearby ultra -luminous x-ray source (ULX) 
XMMU J004243. 6+412519 (iMiddlet.on et al.ll2013h . Sec- 
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tion [5] describes details of the data selection and analysis 
procedure which were employed to produce the results pre- 
sented in section [21 Section U discusses the implications of 
this detection in the context of HE and VHE 7-ray emission 
from BL Lac objects. 



2. Data Analysis 

The Fermi-h AT data were analysed using version v9r27pl 
of the Fermi Science ToolfQ in conjunction with the 
P7S0URCE_V6 (Pass 7) instrument response functions. The 
dataset under consideration was extracted from a circular 
region of interest (ROI) of radius 10 degrees centred on 5C 
3.178 and comprised all SOURCE class photon-like events de- 
tected by the Fermi-LAT before MJD 56282 (2012-12-21) 
with energies between 200 MeV and 300 GeV. Good time 
intervals (GTIs) were generated using the recommended se- 
lection expressiorQ and ROI-based maximum zenith angle 
culEI, before being applied to isolate subsections of the data 
for subsequent analysis. 

The reduced data were analys ed using an un -binned 
maximum likelihood approach (Ma ttox et aLlflQQGl ) imple- 
mented by the gtlike utilitjo, in which two alternative 
model hypotheses were compared by maximising their re- 
spective likelihoods with respect to the the observed photon 
distribution. As a nominal null hypothesis, the expected 
7-ray signal within the ROI was modelled using a combi- 
nation of the standard galactic (gal_2yearp7v6_v0) and 
isotropic (iso_p7v6source) diffuse emission models. The 
model also incorporated all point-like 7-ray emitters within 
20° of 5C 3.178 that are hsted in the 2FGL catalogue. To 
prevent genuine variability or statistical fluctuations of the 
signals from nearby 7-ray emitters from affecting the anal- 
ysis, the overall flux normalisations of sources within a 10° 
radius of the target were treated as variable parameters. In 
order to ameliorate any local mis-modelling of the diffuse 
7-ray emission within the ROI, the overall normalisations 
of the isotropic and galactic template components were also 
allowed to vary during the likelihood analysis. To limit the 
number of degrees of freedom in the final likelihood, the 
remaining source model parameters were fixed to the best 
fitting values published in the 2FGL. A second, alternative 
hypothesis included a model for putative emission from an 
additional point-like 7-ray source with a power-law spectral 
shape parameterised as 



dN/dE = No{E/Eo) 



-r 



where Lmax,o a-nd £max,i are the maximum likelihood val- 
ues obtained when fitting the observed data using the null 
and alternative hypothesis models respectively. If the null 
hypothesis is true, then yjTS is approximately equivalent to 
the source detection significance in the framework of stan- 
dard normal theory. 

3. Results 

As explained in fJU the initial source detection corresponded 
to Fermi data that were obtained during the 200 day flaring 
cycle of the ULX XMMU J004243.6+412519. A subsequent 
analysis of the complete Jermi-LAT dataset revealed strong 
evidence {TS = 54.1 « 7.4o-) for a spectrally hard (F w 1.7) 
point-like HE 7-ray source close to the position of 5C 3.178, 
with a time-averaged 0.2-300 GeV flux of (1.78±0.79) x 10"^ 
ph cm~^s~^. Given its apparent brightness, the absence 
of this source fro m both th e first and second Fermi-hAJ: 
source catalogues (|Abdo et a l. 2010; Nolan et al. 2012) sug- 
gests a transient origin for the observed 7-ray emission. In 
the absence of any a-priori expectation for likely variabil- 
ity timescales, the overall dataset was subdivided into 200 
day intervals (labelled A-H) to match the duration of the 
initial detection dataset. Each interval was individually 
re-analysed to obtain the results listed in Table [T] and illus- 
trated in Figured] 

Assu ming that the test sta tistic is approximately Xi dis- 
tributed (Ma ttox et al.l [19961 ) . then the probability of ob- 
taining a value of TS which exceeds some threshold TSq in 
a single trial is 



Pi{TS > TSq) 



xl{x)dx 



(3) 



TS 



Under the reasonable assumption that the measured 7-ray 
signals corresponding to non-overlapping time slices are sta- 
tistically independent, then the probability of obtaining a 
TS > TSo at least once after testing all 8 sub-intervals is 



Ps 



(4) 



Defining <I>^^ to be the complement of the quantile of the 
Xi distribution, Pg can be used to derive a trials-corrected 
test statistic 



(1) TS'^^-\Ps 



(5) 



The normalisation (A^o) and spectral index (F) were allowed 
to vary during the likelihood optimisation, while the energy 
scale (Eq) was fixed to 1 GeV. 

Evidence for the detection of an additional HE 7-ray 
source was evaluated in terms of a likelihood ratio test 
statistic 



A trials-corrected test statistic value of TS' — 58.9 
(~ 7.7cr) provides clear evidence for transient 7-ray emis- 
sion during interval G, which begins on MJD 55882 (2011- 
09-17). Further subdivision of interval G into 10 week, 2 
week and 1 week time slices does not reveal significant fiux 
variability on shorter timescales. The source was not con- 
clusively detected during any other interval, although there 
is marginal evidence for emission (TS' ~ 8.5) during in- 
terval H. For intervals yielding TS < 9, Table [1] lists 95% 
confidence level upper limits on the 0.2-300 GeV fiux, which 

1 http://fermi.gsf c.nasa.gov/ssc/data/ana lysis/scit^^l iTd^n^^m^^ Figure[Il The time-resolved analysis also 

reveals that the source would not have been detectable dur- 
ing the first 24 months of the Fermi mission, which explains 
its absence from the Fermi-LAT source catalogues. 

Given the limited duration of the detectable 7-ray emis- 



TS= -2 In 



L 



max.O 



(2) 



(LAT_C0NFIG==1) 



^ (DATA_QUAL==1) 
ABS(R0CK_ANGLE)<52 

^ This cut excludes all time intervals when the zenith angle of 
any part of the ROI exceeds 100° 



|http: //f ermi.gsf c .nasa.gov/ssc/data/analysis/scitools/Isaibp/gSilitfeeqtt emt analyses only considered the restricted 
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Interval 


^start 

(MET) 


^cnd 

(MET) 


^start 

(UTC) 


TS 


TS' 


Flux (0.2-300 GeV) 
{10-^ ph cm-^s-i) 


r 


A 


239557417 


256837417 


2008-08-04 15:43:36 


0.0 


0.0 


< 2.39 




B 


256837417 


274117417 


2009-02-20 15:43:35 


3.0 


0.4 


< 5.37 




C 


274117417 


291397417 


2009-09-08 15:43:35 


3.3 


0.6 


< 5.80 




D 


291397417 


308677417 


2010-03-27 15:43:35 


4.1 


1.0 


< 1.36 




E 


308677417 


325957417 


2010-10-13 15:43:35 


6.8 


3.2 


< 7.81 




F 


325957417 


343237417 


2011-05-01 15:43:35 


7.1 


3.5 


< 4.46 




G 


343237417 


360517417 


2011-11-17 15:43:35 


63.0 


58.9 


8.22 ±2.04 


1.76 ±0.09 


H 


360517417 


377797417 


2012-06-04 15:43:35 


12.3 


8.5 


3.05 ± 1.32 


1.69 ±0.12 


All 


239557417 


377797417 


2008-08-04 15:43:36 


54.1 




1.78 ±0.79 


1.70 ±0.15 



Table 1. Analysis results corresponding to eight 200 day subintervals (A-H) as well as the complete Fermi-LAT dataset. Test 
statistic values are quoted before (TS) and after {TS') correcting for the number of independent trials. For intervals during which 
the source was detected such that TS > 9, the HE 7-ray flux and best fitting power-law spectral index (F) are also listed. For 
other intervals, upper limits to the HE 7-ray flux (at the 95% confidence level) are presented. 
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Fig. 1. Left panel: Values of TS corresponding to intervals A-H. Right panel:. 95% confidence upper limits on the 0.2-300 GeV 
7-ray fiux are presented for intervals during the intervals in which the source was not detect ed. Th e shaded areas in both plots 
indicate the temporal extent of the Fermi-LAT dataset used to compile the 2FGL catalogue l|Nolan et al.,,2012i ). 



dataset of events detected during interval G, in order to 
maximise the signal-to-noise ratio. The gtf indsrc utilitj(f| 
was applied to derive the most likely coordinates of the ob- 
served 7-ray emission, yielding (aLAT,j20oo, ^lat,J20oo) = 
(00M8'"20''.89,+39°52'26".l) with a 95% confidence er- 
ror circle radius of 0.06°. Using these coordinates, 
the automatic so urce association method described by 
(jAbdo et al.l |2010( ) and implemented by the ScienceTools 
executable gtsrcid, was applied to search for plausible 
multi-wavelength counterparts to the observed HE 7-ray 
emitter. After considering all sources listed within the cat- 
alogues specified in Table 9 of iNolan et al.l ([2012 '). the pro- 
cedure identified 5C 3.178 as the only viable candidate, de- 
riving a posterior probability of 91% that the association is 
correct. 



The Fermi-hAT data corresponding to subinterval G 
were used in conjunction with the gttsmap too^l to gener- 
ate the TS sky map presented in the left panel of Figure 
[5] To generate this map, the null hypothesis model de- 
scribed in ^J2]was used, with all parameters except for the 
diffuse component normalisations fixed to the best fitting 
values that were obtained for this subinterval during the 
likelihood analysis. The transient 7-ray source is clearly 
visible as a point-like enhancement of the computed TS 
values, coin cident with the nominal position of 5C 3.178 
(|Maoll2011h . 

The observed 0.2-300 GeV spectral energy distribution 
and best-fitting power law model corresponding to interval 
G are presented in the right-hand panel of Figure [5J The 
data correspond well with the assumed model and indicate 



|http: //fermi .gsfc ■ nasa.gov/ssc/dat a/ analysis/scitools/hfe!|]pl^-gif/!l^flgnTTTri 

Article number, page 3 of [5] 




Fig. 2. Left panel: Unbinned likelihood test statistic map corresponding to subinterval G of the Fermi-ljAT dataset. The white 
cross marks the nominal radio position of 5C 3.178. Right panel: HE 7-ray spectral energy distribution between 0.2 and 300 GeV 
corresponding to subinterval G of the Fermi-LAT dataset. The dotted line indicates the best-fitting power-law model, and the 
grey region indicates the la fit uncertainties. Residuals of the observed data points with respect to the best fitting power law 
spectral model are also presented. 



that the transient 7-ray emission is spectrahy hard, with a 
power law index To = 1.76 ± 0.09. 

4. Discussion 

The results presented in SJ2] provide clear evidence for tran- 
sient, spectrally hard HE emission from a point-like 7-ray 
source that is positionally coincident with the BL Lac ob- 
ject 5C 3.178. The spectral properties and temporal vari- 
ability of the transient are both consistent with the char- 
acteristics of other BL Lac objec ts that have previously 
been detected by th e Fermi-LAT (jAckermann et ahl 120111 : 
l^entiirk et al■l[2013^ . On the basis of this combined evi- 
dence, it seems reasonable to attribute the observed 7-ray 
emission to a HE flare of the BL Lac object 5C 3.178. 

The hard spectral index of the observed HE emission 
is also consistent with the sub-population of 36 High- 
Synchrotron-Peaked BL Lac objects (HBLs) that have also 
been det ected at VHE 7-ray energies, all of which have 
F < 2.3 (Igentiirk et al.ll2013f ). Indeed, the moderate red- 
shift of 5C 3.178 suggests that the source may be a vi- 
able target for observations by ground-based Cherenkov 
telescope arrays. Combined HE and VHE 7-ray obser- 
vations HBLs have been used to probe the p roperties of 
the extragalactic background hght (e.g. lAckermann et al.l 
l2012t IH.E.S.S. Collaboration etaIll2013^ and can provide 
strong constraints on the value of bulk cosmological prop- 
erties s uch as the strength of the extrag alactic magnetic 
field (T aylor et al,, 2011: V ovk et al.ir2012l ). 

Studies of this type invariably benefit from multi- 
wavelength observations and strictly simultaneous data are 



preferred when dealing with variable and unpredictable 7- 
ray signals. Without an a-priori expectation of transient 
7-ray emission, achieving this simultaneity requires rapid 
identification of potential targets in order to coordinate the 
instruments involved. Atmospheric Cherenkov telescope ar- 
rays such as MAGI C II (Aleksic et al. 2012) or VERITAS 
(jHolder et al.l[2006h provide good sensitivity above ~ 100 
GeV, but are characterised by limited fields-of-view (< 5°^) 
which renders the serendipitous discovery of flaring 7-ray 
sources using these instruments unlikely. In contrast, con- 
tinuous survey observation mode of the Fermi-hAT effec- 
tively guarantees that any transient emission lasting longer 
than ~ 6 hours will at least be observed, although limited 
sensitivity in the GeV energy range may lead to fainter 
sources escaping detection. 

Systematic blind searches for variable signals within the 
Fermi-LAT all-sky dataset are extremely computationally 
intensive and can only consider a limited number prese- 
lected variability timescales. Currently, automatic moni- 
toring for transient 7-ray signals in the latest Fermi-LAT 
data is provided by the Fermi Flare Advocate (FA) ser- 
vice (|Ciprini fc Fermi-LAT Collaboration! I2012D . This fa- 
cility performs preliminary analyses that are optimised to 
search for bright {Fe>womcV ^ 10~^ phs^^'^cm^) sporadic 
7-ray emission within intervals lasting 6 hours, one day, 
and one week. Restricting transient searches in this manner 
inevitably sacrifices sensitivity to signals with longer vari- 
ability timescales. Indeed, if the proposed association with 
5C 3.178 is correct the 7-ray signal identified in this work 
may indicate a population of mildly variable BL Lac ob- 
jects that are overlooked by the current FA policy. Assum- 
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ing that this is the case, the temporal characteristics of the 
observed HE signal appear to motivate all-sky searches for 
moderately bright 7-ray emission that varies on ~ 100 — 200 
day timescales. 

5. Summary 

Fermi-LAT data have been used to identify a point-like 
source of transient HE 7-ray emission that is positionally 
coincident with the BL Lac object 5C 3.178. The probabil- 
ity that 5C 3.178 is the correct association for the transient 
emission was calculated to be 91%, while a search for alter- 
native counterparts yielded no viable candidates. A tempo- 
rally resolved analysis yielded a trials-corrected maximum 
likelihood ratio test statistic of 58.9 (equivalent to sa 7.7it) 
and localised the observed 7-ray signal to the time interval 
between MJD 55882 (2011-11-17) and MJD 56102 (2012- 
06-24). During this interval, the source was characterised 
by a hard power law energy spectrum (dN/dE oc E~^) 
with spectral index T = 1.76 ±0.09 and a 0.2-300 GeV flux 
-Fb.2-300GcV = (8.22 ±2.04) x 10"^ ph cm-^g-i. In combi- 
nation, these spectral and temporal properties are reminis- 
cent of a limited number of HBLs that have previously been 
observed by the Fermi-L AT , and which are often VHE 7- 
ray emitters. The detection of 5C 3.178 using ground-based 
atmospheric Cherenkov telescope arrays would extend the 
existing catalogue of BL Lacs that have been detected at 
both GeV and TeV energies, enabling a more secure clas- 
sification of the source and simultaneously establishing 5C 
3.178 (z « 0.25) as the fourth most distant object in this 
categorjQ. As well as elucidating the intrinsic processes 
that operate within AGNs, simultaneous HE and VHE ob- 
servations of distant HBLs can also be used to constrain the 
properties of extragalactic photon fields or the intergalactic 
magnetic field. 
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